Draft version January 15,2013 

Preprint typeset using l^TgX style emulateapj v. 1 1/10/09 



ASYMMETRY AND THE NUCLEOSYNTHETIC SIGNATURE OF NEARLY EDGE-LIT DETONATION IN WHITE 

DWARF CORES 

David A. Chamulak'-^, Casey A. Meakin^--', Ivo R. Seitenzahl'*, James W. Truran'-^-^-^ 

Draft version January 15, 2013 



o 

(N 
D 



00 
O 



> 

O 

rn 

rn 
O 



ABSTRACT 

Most of the leading explosion scenarios for Type la supemovae involve the nuclear incineration of a white 
dwarf star through a detonation wave. Several scenarios have been proposed as to how this detonation may 
actually occur, but the exact mechanism and environment in which it takes place remain unknown. We explore 
the effects of an off-center initiated detonation on the spatial distribution of the nucleosynthetic yield prod- 
ucts in a toy model - a pre-expanded near Chandrasekhar-mass white dwarf. We find that a single-point near 
edge-lit detonation results in asymmetries in the density and thermal profiles, notably the expansion timescale, 
throughout the supernova ejecta. We demonstrate that this asymmetry of the thermodynamic trajectories should 
be common to off-center detonations where a small amount of the star is burned prior to detonation. The sensi- 
tivity of the yields on the expansion timescale results in an asymmetric distribution of the elements synthesized 
as reaction products. We tabulate the shift in the center of mass of the various elements produced in our 
model supernova and find an odd-even pattern for elements past silicon. Our calculations show that off-center 
single-point detonations in carbon-oxygen white dwarfs are marked by significant composition asymmetries in 
their remnants which bear potentially observable signatures in both velocity and coordinate space, including 
an elemental nickel mass fraction which varies by a factor of two to three from one side of the remnant to the 
other 

Subject headings: nuclear reactions, nucleosynthesis, abundances — supernovae: general — white dwarfs 



1. INTRODUCTION 

The fact that their peak absolute magnitude is correlated 
with the width of the light curve has allowed Type la super- 
novae (SNe la) to be used as standard candles in determin- 
ing cosmological distances. Despite this widespread use of 
SNe la as standard candles, many problems in describing how 
the explosion of the star happens still remain. 

A number of models ha ve been proposed to explain 
SNe la (for a review, see iHiUebrandt & N iemevei^ 120001 
[Ppdsiadlowski et al. 200^. Two promising candidates for 
an explosion mechani sm are the sub-Chandras e khar dou- 
ble d etonation model (i Wooslev & Weave3 Il994t iFink et al] 
120101) and the gra vitational confined detonation (GCD) model 
dPlewa etaT] 1200 4). In both models, most if not all of the 
nuclear burning occurs in a detonation wave. In the sub- 
Chandrasekhar model a layer of helium is deposited on the 
surface of a white dwarf. The helium layer detonates re- 
sulting in a shock wave traveling around the surface of the 
white dwarf. When the shock wave converges at the an- 
tipode, a detonation is t hought to be triggered off center in 
the carbon-oxygen core (IWooslev & Weaverll 19941 iFink et all 
[TOIO; Sim et al. 20ldi). The nucleosynthetic yield is set pre- 
dominantly by the mass of the carbon- oxygen white dwarf 
core and the mass of the helium layer (IWoosley & Weaver! 
n^lFiiiketaniMol) . 
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In the GCD model, the carbon burning runaway within the 
convective core of a near Chandrasekhar-mass white dwarf is 
postulated to occur in a small region displaced from the stel- 
lar center that becomes a highly buoyant flame bubble and 
quickly rises to the stellar surface after b urning only a few 
percent of the star during its ascent (e.g. iPlewa et aP 120041 
iLivne et al. 2005). When the buoyant ash, as well as unburned 
material pushed ahead of the rising flame bubble, erupts forth 
from the stellar core it is largely confined to the surface of the 
white dwarf by gravity. It then becomes a strong surface flow 
that sweeps completely over the star, eventually converging 
in a region opposite to the breakout location. Although the 
details of the detonation initi ation process i s still a topic of 
active research (see e.g., Ro pke et al.ll2007t ISeitenzahl et alJ 
2009a,b) the high temperatures and densities reached within 
the converging surface flow are strong indicators that a deto- 
nation is likely to be triggered. 

The resulting GCD nucleosynthetic yield consists almost 
entirely of detonation burning products and depends on how 
much the star has expanded by the time the detonation ini- 
tiates. More highly expanded (hence lower density) cores at 
detonation result in a smaller fraction of Fe peak nuclei, less 
^^Ni, and consequently a larger fraction of intermediate mass 
elements (IMEs) due to incomplete relaxation to nuclear sta- 
tistical equilibrium (NSE). Therefore, lower luminosity (less 
^^Ni producing) explosions are accompa nied by a large r yield 
of IMEs, as has been observed (see e.g., Mazzali et al.f 2007). 

The expansion of the star prior to detonation, in a GCD 
scenario, results from the work done by the rising flame bub- 
ble. Deflagrations that bum more mass prior to reaching the 
stellar surface excite higher amplitude pulsations and hence 
more-expanded stars at the time of detonation. It has been 
found that the expansion of the star due to the deflagration 
is very well represented by the fundamental radial pulsation 
mode of the underlying white dwarf (see Fig. 2 and Fig. 12 in 
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iMeakin et alj|2009h . Therefore, while it is essential to under- 
stand the nature of the deflagration so as to better understand 
the mapping between initial conditions and final outcomes, 
the range of outcomes due to the deflagration can be parame- 
terized by the pulsation amplitude, resulting in a one parame- 
ter family of models. While the phase of the pulsation at the 
time of detonation is an additional parameter, it plays a lesser 
role since the pulsation period is longer then the detonation 
crossing time. 

In this Paper we explore the nucleosynthetic yields that 
result from an edge-lit detonation of a pre-expanded near 
Chandrasekhar-mass white dwarf core out of hydrostatic equi- 
librium. For our purposes we define edge-lit detonation to 
be an off-center detonation where very little of the star has 
burned before hand. This is a toy-model which neglects the 
effect the deflagration has on the nucleosynthesis. However, 
due to the small amount of mass burned in the deflagration it 
captures most of the essential features of the GCD model. The 
physical mechanism that gives rise to the abundance asym- 
metries in the ejecta will furthermore manifest itself in other 
off-center detonation models, such as the sub-Chandrasekhar 
double detonation models. 

In Section|2]we describe the explosion model and our com- 
putational method. Section [5] is comprised of three parts. 
First, we discuss how the dynamics of an asymmetric detona- 
tion affect the hydrodynamic profile of the SN la remnant. We 
then consider how variations in the thermodynamic trajecto- 
ries of detonated material affect the resulting nucleosynthesis. 
Finally, we present the elemental yields found for our explo- 
sion model and quantify the asymmetric distribution thereof 
in terms of center of mass offsets, which we provide in tabu- 
lated form. We also show that there is an asymmetric distribu- 
tion in velocity space for elemental nickel. We then conclude 
in Section HI with a discussion of how our results relate to 
previous work and how our work can be improved. 

2. THE EXPLOSION MODEL 

The explosion model discussed in this paper involves the 
detonation of a cold (T=4 x 10^ K) white dwarf of mass 
1.365Me, comprised of 50% '^C and 50% '^O by mass, 
which has been expanded according to its fundamental radial- 
pulsation mode by such an amplitude that it has a central den- 
sity of 10** gcm^^. A detonation is initiated in this expanded 
white dwarf at a radial location of r = 2 x lO^cm, where the 
density of the white dwarf is p = lO^gcm^^, by heating a 
small spherical volume with radius rini ^ 4 km to high tem- 
perature Tini ^ lO^K. This heating immediately triggered a 
detonation which propagates away from the point of initia- 
tion. 

The subsequent reactive-hydrodynamic evolution of the 
detonated white dw arf was conducted using the FLASH code 
(iFryxell et a l. 2000*). The code frame work and the include d 
physics is identical to that described in iMeakin et al.l (l2009h . 
The star was set up on a Cartesian grid in 2D with cylindri- 
cal symmetry. An effective adaptive mesh refinement resolu- 
tion of 1 km is used. The detonation was initiated along the 
symmetry axis, which is the only natural way to capture its 
evolution in this geometry. Due to the cylindrical symmetry 
of the physical scenario and the lack of large scale instabil- 
ity in the detonation front this setup gives the s ame results as 
a full 3D calculation of the same problem (see IMeakin et al.l 
( 120091) which provides an in depth description of the detona- 
tion phase, and a comparison between 2D and 3D). 
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Si 
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Ar 



17-28 
20-31 
20-33 
22-35 
22-38 
26-40 
27^2 
31^4 
31^7 



K 

Ca 
Sc 
Ti 
V 
Cr 
Mn 
Fe 
Co 



35^6 
35-53 
40-53 
39-55 
43-57 
43-60 
46-63 
46-66 
50-67 



Ni 
Cu 
Zn 
Ga 
Ge 
As 
Se 
Br 
Kr 



50-73 

54- 70 

55- 72 

58- 73 

59- 76 
62-76 
62-82 
71-81 
71-86 



Energy release and bulk compositional evolution due to nu- 
clear burning is coupled to the hydrodynamic flow by ad- 
vancing a system of progress variables which represent three 
stages of burning: (1) the burning of carbon and oxygen to 
silicon group elements, (2) the relaxation of the silicon group 
elements to NSE, and (3) the evolution of material which has 
already relaxed entirely to NSE, including its adjustment to 
changing density and temperature and neutronization due to 
weak interactions among the Fe peak elements. Additional 
details can be found in ICalder et al.l (l2007l) : IMeakin etaP 
(.2009). 

The pre-expansion of the white dwarf prior to detonation 
initiation results in IMq of high density material (p > 
10^ gcm^^), which burns to NSE in the detonation (primar- 
ily as 56Ni). The remaining ^ O.SlMc., of material bums to 
IMEs, e.g. Si, S, Ca, and Ar, resulting in only a small amount 
(< O.O4M0) of unbumed C/O in the outermost layers of the 
remnant. 

Detailed yields are calculated by post processing La- 
grangian tracer particles included in the explosion calcula- 
tion and are the primary focus of this paper We recorded the 
time history of 10^ particles which were initialized to evenly 
sample the initial mass of the white dwarf. Our reaction net- 
work incorporates 493 nuclides from n to ^^Kr (Table[T]i. We 
use the reaction rates from the J oint Institute for Nu clear As- 
trophysics REACLIB Databas4!l (IC yburt et ani201(X and ref- 
erences therein); the light-element rates are mostly experi- 
mental and are from compilations such as'Caughlan & Fowlerj 
(1988) and Iliadis et al. (2001). Weak reaction rates are taken 
from [ Fuller et al.l (I1982D and Lan ganke & Martmez-Pinedol 
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(120011). Screening is incorporated using the formaUsm of 
IGraboskeetaHdigTl . 

3. RESULTS 

3.1. Explosion Dynamics and Remnant Asymmetry 

The detonation propagates from the point of initiation 
at nearly the Chapman-Jouguet (CJ) speed, Dcj ^ 1.2 x 
10^ cms"' at a Mach number of Mcj ^ Dcj/c, ^ 3.4. Be- 
cause of the weak upstream density dependence of the deto- 
nation speed under these conditions, the detonation front re- 
mains very nearly spherical in shape as it engulfs the star 
The total time required for its passage across the expanded 
white dwarf core is tcmss ^ '^f/Dcj ^ 0.5 s. This is followed 
by a period of ^ 0.5 s in which the pressure forces drive 
the completely incinerated remnant into a homologous expan- 
sion, characterized by a purely radial expansion velocity pro- 
file with an expansion rate proportional to the radial position 

^ |http : //groups ■ nscl .msu ■ edu/ j Ina/reacllb/db/ | 
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V °^ r. After only a few seconds, the remnant is expanding 
ballistically, and the total energy budget is dominated by the 
kinetic energy. The homologous velocity profile results in a 
self-similar density profile which persists until the remnant 
begins to interact with the interstellar medium. 

The expanding remnant resulting from the detonation is 
marked by significant asymmetry. The late time density pro- 
files along the symmetry axis and the equator are shown in 
Fig [T] scaled by the peak density in the remnant pc- The 
initial, spherically symmetric white dwarf density profile is 
also shown for comparison. It can be seen that the density 
peak is shifted into the hemisphere in which the detonation 
was initiated, y > in this case, resulting in a steeper density 
gradient in this hemisphere. This is i n agreement with th e se- 
ries of GCD simulations described in lMeakin et al.l(l2009 ')(see 
their Fig. 9-11) which show that the density isocontours 
in the remnant are well described by concentric circles that 
have centers offset from the initial stellar center The density 
isocontours were found to have larger offsets at higher den- 
sities, with the largest offset centered on the peak density in 
the remnant. In all cases the density peaks on the side of the 
remnant where the detonation originated and has a steeper 
density gradient in that region. 




0.01 i I . . . .: 

-2-10 1 2 

scaled offset position 




-2-10 1 
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Fig. 1. — Late time (t > 3 s) density and velocity profiles for the post- 
detonation state. The density is scaled by the peak value and the position is 
scaled by the density e-folding distance in the equatorial direction. The thick 
gray line shows the scaled density profile of the initial white dwarf" model, 
while the post-detonation state model is shown by thick black lines for lines 
through the equator (solid) and through the poles (dashed). 

UnUke the density profile, the velocity profile (also shown 
in Fig.[T]) does not show an asymmetry, but is everywhere ra- 
dially directed and spherically symmetric. This leads to an 
asymmetry in the density as a function of expansion veloc- 
ity, which is likely to result in a viewing angle dependence 
for the Ught curve and the spectral signature. Related com- 



position asymmetries, discussed in ii i]3.2l and l3.3l below. also 
contribute to observable asymmetries and viewing angle de- 
pendencies. 

A revealing format for presenting the dynamics of the det- 
onation and the subsequent expansion is the space-time dia- 
gram. In Fig.|2]we present the space-time trajectories for all 
of the tracer particles that were initialized near the symme- 
try axis of the white dwarf. The left panel shows the time 
period over which the detonation traverses the stellar core, 
while in the right panel we show the later time evolution that 
ends in a radially expanding, ballistic trajectory for each of 
the particles. The bold dashed line in the left panel shows the 
path taken by a theoretical detonation having a constant speed, 
which matches the kinks in the particle trajectories very well. 

The following features in this figure merit further discus- 
sion. (1) The trajectories are slowly converging prior to det- 
onation. This is the signature of the stellar core undergoing 
mild contraction as a result of having been expanded by a ra- 
dial pulsation mode prior to detonation. (2) The detonation 
accelerates material in the direction it is propagating. This is 
the primary source of the asymmetry imprinted on the rem- 
nant at late times. A large number of Lagrangian tracer parti- 
cles in the detonated hemisphere are first accelerated towards 
the stellar center by the detonation before they are turned 
around by pressure forces and accelerated to their final, out- 
wardly directed expansion velocities. The exact number of 
tracer particles accelerated towards the stellar center by the 
detonation is dependent on distance from the center of the star 
to the point where the detonation was initiated, a. A given 
tracer particle with a central angle and distance from the 
center of the star r will be accelerated towards the center of 
the star if ^^(Q) >r/a. On the other hand, tracers in the oppo- 
site hemisphere are accelerated by the detonation in the same 
direction as their final expansion velocity, reaching their final 
velocity on a shorter timescale. (3) The material in the det- 
onated hemisphere is accelerated to lower velocities overall 
compared to material in the opposite hemisphere (see Fig. |2] 
(right)). This mapping between initial position (and therefore 
initial density) and resultant expansion velocity explains the 
density profile: the material lines in the more rapidly expand- 
ing hemisphere are stretched out over a larger region of space, 
and hence to a lower relative density, than the more slowly 
expanding regions. (4) A natural consequence of the explo- 
sion dynamics is an asymmetry in the expansion timescale 
texp, defined as the time required for the detonated material 
to drop from its post-detonation temperature maximum T,nax 
to e^^Tmax, resulting from the differential rate at which ma- 
terial cools (nearly adiabatically) due to the post-detonation 
expansion. This follows directly from point (2) above. The 
expansion timescale asymmetry is shown in Fig. [3] where we 
have plotted the tracer particles at their initial position in the 
stellar core, color coded by their post-detonation expansion 
timescale. It is obvious from this figure that the material in the 
detonated hemisphere [y > 0) has overall a larger expansion 
timescale than in the opposite hemisphere f or a give n in itial 
upstream density. As will be discussed in i^ il3.2l and l3.3l this 
expansion timescale distribution imparts an asymmetry in the 
resultant nucleosynthetic yield. 

3.2. Nucleosynthesis Dependence 

We find that nuclear burning in SNe la progresses in three 
distinct stages (iKhokhlovll 1 983L 1 1 99 11) . The first stage is cai-- 
bon burning. During carbon burning '^C-i-'^C is the primary 
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Fig. 2. — (left) On-axis Lagrangian tracer particle positions are shown as a function of time in this space-time diagram. The thick dashed line shows the 
theoretical position of a constant speed detonation originating from the location y = 0.2 x lO' cm, coincident with the detonation initiated in the simulation. A 
theoretical detonation speed taken to be I'j^t = 11 x lO' cms"' for the detonation as it moves inward toward the high density core and a speed that is 10% lower 
as it moves outward into the low density surface material provides a very good match to the simulation data, (right) At times greater than Is the Lagrangian tracer 
particles exhibit homologous expansion. 




2 4 6 8 

X (10^ cm) 

Fig. 3. — Initial spatial position of the Lagrangian tracer particles. Color 
represents each particle's expansion time scale. The center of the star is at 
X=0 Y=0 The detonation was initiated at X=0 Y=2. Particles on the side of 
the star where the detonation starts have higher expansion time scales than 
the particles on the opposite side of the star. The expansion time scale is 
calculated from the temperature profile of each tracer. 



reaction taking place. We find that carbon burning never 
reaches an equiUbrium state in a small region of the star where 
the final carbon mass fraction is above 10^"*. The carbon burn- 
ing reactions are sensitive to temperature with 

^^^c</(r9)r-2/V«4-i^<''. (1) 

dt 

Here dY {^^C) / dt is change in '^C abundance over change in 
time, Tg is temperature in units 10^ K, and f{Tg) is a func- 
tion defining the temperature effect on the branching ratio 
between th e '^C(^^C , a)^°Ne and ^^C(}^C, pf^Na reaction 
(|Caughlan & FowleJ fl988). Any change in the thermal pro- 
file will result in a different abundance pattern for material 
that does not proceed to the next phase of burning. 

The next stage is oxygen burning. Here '^O and the prod- 
ucts of carbon burning proceed to silicon group elements, like 
S, Ar, and Si. As in carbon burning, oxygen burning never 
reaches an equilibrium state and therefore also shows a de- 
pendence on the thermal history. We find X('^O) > 10^^ if 
the next burning stage did not start. Very little mass of the 
star (< O.O4M0) is in a region that does not complete either 
carbon or oxygen burning, so most of the star proceeds to the 
next burning stage. 

At higher temperatures and densities, silicon burning is the 
dominant form of nucleosynthesis. In silicon burning groups 
of nuclides enter into equilibriu m; a state known as quasi- 
statistical equilibrium (QSE; e.g.. lWooslev et al.ll973h . There 
are two ways in which changing the expansion time scale can 
affect the abundances in QSE. First, while equilibrium holds 
within a group of nuclides, it does not hold between groups. 
Second, within a group of nuclides in equilibrium reactions 
freeze out at different temperatures, resulting in an abundance 
pattern that depends on the expansion time scale. This will be 
further explored below. 
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Fig.|4]shows how the abundances of silicon and nickel vary 
as functions of expansion time scale over a small range of 
temperatures for tracer particles that never finish silicon burn- 
ing. Even with the scatter from plotting particles with differ- 
ent peak temperatures, there is a clear dependence of the abun- 
dances on expansion timescale. The directions of these trends 
are counterintuitive, but Fig. |5]shows their origin. Both tracer 
particles shown reach a peak temperature of 5 x 10^ K, and 
their nucleosynthesis is nearly identical up to that point. The 
particle with the longer (0.6016 s) thermal expansion time 
scale, however, has a temperature that falls faster over the first 
0.1 s than the particle with the shorter (0.4428 s) thermal ex- 
pansion time scale. This rapid decrease in temperature results 
in less ^^Si burned to ^^Ni even though the thermal expansion 
time scale is longer. 




0.0 1 , , , \ , . , \ , . , I 

0.2 0.4 0.6 0.8 

expansion timescale (s) 

Fig. 4. — Mass fraction of ^^Si(red) and '^Ni(blue) as a function of expan- 
sion time scale. This plot was made from tracer particles that had a maximum 
temperature between 4.99 x lO' K and 5.01 x lO' K and a density of approx- 
imately 1.5 X lO' gem"'. The lines shown are least squares fits to the data. 

Material exposed to high enough density and temperature 
conditions for long enough will arrive in a state of nuclear sta- 
tistical equilibrium (NSE). In this state all nuclear reactions 
enter equilibrium and lose all history of the thermal evolu- 
tion up to that point. For material that has reached this state, 
the only dependence that the final yield has on expansion 
timescale occurs during the process of freeze out. Freeze out 
occurs for a nuclide when the temperature drops low enough 
that all strong reactions become too slow to change the nu- 
clides abundance again. Because this condition occurs at dif- 
ferent temperatures for different nuclides, the final yield de- 
pends on the rate at which material cools. Like in QSE, dif- 
ferent reactions freeze out at different temperatures, leading 
to yields that depend on the expansion timescale. Therefore, 
all three stages of burning and NSE in SNe la are affected 
by different expansion time scales, resulting in a clear com- 
positional asymmetry that will be discussed in the following 
section. 

3.3. Phenomenology 

We now present the results of our reaction network calcu- 
lations for a near edge-lit detonation in a SN la model. We 
find that a number of nuclides exhibit pronounced asymme- 
tries across the stellar remnant. We quantify this effect by 
calculating the center of mass for a given element. Suppose 



a tracer particle / has a mass fraction X,(Z) of element Z, and 
its position is given by vector r,. Then the center of mass for 
a given element, rc„,(Z), is given by the equation 

r,,„(Z) ^ (^X,iZ)r}j I {^X,{Z)]^ . (2) 

Due to the cylindrical symmetry of our explosion model the 
displacement of the center of mass for any element lies along 
the y axis. A velocity for the center of mass v„n can be calcu- 
lated by replacing r,- by the velocity v; in equation|2l Table|2] 
shows the total mass, displacement of the center of mass, and 
velocity of the center of mass for elements between carbon 
and germanium. For reference, the detonation was initiated at 
^ 2 X 10** cm. These numbers are correct for the end of our 
simulation (f ^ 3 s), where strong reactions have frozen out 
and homologous expansion has been reached. Some of the 
isotopes making these elements, ^^Ni for example, decay so 
Table |2] evolves with time. The last three columns show the 
total mass, displacement of the center of mass, and velocity of 
the center of mass assuming all radioactive elements instantly 
decayed to their stable isotopes. Note that nickel shows no 
change in center of mass or velocity since most of the mass of 
the star ends up as radioactive nickel in our model. If the com- 
plete star was burned to nickel then by definition the change 
in center of mass would be zero since our model conserves 
mass and momentum. Elements lighter than silicon have their 
masses distributed more in the direction where the detonation 
was initiated. Elements heavier than silicon are, for the most 
part, distributed away from where the detonation was initi- 
ated. These elements also display an odd-even pattern where 
odd Z nuclei, like cobalt and copper, are predominantly dis- 
tributed farther away from the start of the detonation than their 
even Z counterparts, like iron and zinc. This asymmetry is 
due to the different thermal histories of the two sides of the 
star affecting the nucleosynthesis as outlined in ^ ^3.2\ 

Our simulation was not run sufficiently past freeze out to 
allow all beta unstable nuclei to decay. However, if we decay 
the unstable nuclei and group them in elemental abundances, 
we find elemental nickel to have a clear gradient over the star. 
In material not burned to NSE we find the mass fraction of el- 
emental nickel between +90° and —90° central angle to each 
tracer particle to increase by a factor of 2-3. Figure 0shows 
the dependence of elemental nickel on central angle to each 
tracer particle. 

An interesting side effect of the different expansion times is 
that material on opposite sides of the star expands at different 
velocities. This leads to a gradient in velocity space. Figure|2] 
shows how elemental nickel, iron, manganese, and chromium 
vary with radial velocity for different central angles. In mate- 
rial not burned to NSE, the part of the remnant with the most 
nickel is also the part with the highest radial velocity. This is 
self-consistent since the nickel mass fraction is greater on the 
side of the remnant with the shortest expansion time. There- 
fore, it follows it should have the highest radial velocity. 

4. DISCUSSION 

We have computed the abundances and spatial distri- 
butions of nuclides in an explosion of an expanded near 
Chandrasekhar-mass white dwarf resulting from an off-center 
initiated detonation, a toy model that captures the thermody- 
namic profile of some SN la explosion models. We find a 
compositional asymmetry in the ejecta produced by the det- 
onation. This compositional asymmetry is connected with 
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Fig. 5. — (left) Temperature profile for two particles with similar peak temperatures but different expansion time scales. The expansion time scale for the solid 
line is 0.4428 s and the expansion time scale for the dashed line is 0.6016 s. The times for both particles have been offset such that the peak temperature is 
reached at 1 x 10"^ s. (right) The mass fraction of '^C, ~'*Si, and ^^'Ni as a function of time for the same time-adjusted particles. 



TABLE 2 

Centers of mass and velocities for various elements 





Ac„,(~3 s) 
(lO** cm) 


K-„,(~3s) 
(10'*cms"') 


mass(~ 3 s) 
(g) 


A„„ {decayed) 
(lO'* cm) 


Vci„{decayed) 
(10**cms"') 


msiss{decayed) 
(g) ■ 


c 


29.1 


9.9 


9.00x10^'' 


29.1 


9.9 


9.00x1029 


N 


18.2 


6.2 


2.91x10^5 


11.6 


4.1 


8.41x1026 





7.30 


2.4 


7.28x10^' 


7.30 


2.4 


7.28x10" 


F 


27.9 


9.4 


2.19x10^2 


35.2 


12. 


3.34x1023 


Ne 


28.7 


9.9 


1.03x10'" 


28.7 


9.9 


1.03x1030 


Na 


25.6 


8.8 


1.12x10^^ 


22.3 


7.6 


5.69x102'' 


Mg 


13.8 


4.9 


1.37x10'' 


13.8 


4.9 


1.37xl03' 


Al 


21.6 


7.5 


1.74x10^* 


18.5 


6.4 


3.13x1021* 


Si 


-1.09 


-0.6 


2.18x10'^ 


-1.09 


-0.6 


2.18x1032 


P 


-4.44 


-1.9 


8.34xl02« 


3.77 


1.2 


4.87x1028 


S 


-2.68 


-1.2 


1.18x10'^ 


-2.68 


-1.2 


1.18x1032 


CI 


-9.31 


-3.7 


1.82x102* 


-4.28 


-1.8 


1.21x1028 


Ar 


-2.95 


-1.3 


2.60x10" 


-2.94 


-1.3 


2.60xl03' 


K 


-4.53 


-1.7 


5.12xKF 


-4.20 


-1.7 


3.66xl02^ 


Ca 


-2.55 


-1.1 


2.80x10" 


-2.56 


-1.1 


2.80x1031 


Sc 


-3.36 


-1.3 


2.02x1026 


-9.21 


-3.7 


1.57x1025 


Ti 


-2.78 


-1.1 


3.93 xlO^* 


-2.35 


-1.0 


7.54 xl02' 


V 


-2.04 


-0.7 


3.68 xlO^' 


-3.91 


-1.5 


2.45x1028 


Cr 


-2.37 


-1.0 


7.60x102'' 


-1.81 


-0.8 


1.71x1031 


Mn 


-3.68 


-1.4 


2.58xl02« 


-6.01 


-2.4 


2.91x102' 


Fe 


-1.80 


-0.8 


1.71x10" 


0.00 


0.0 


2.16x1033 


Co 


-5.59 


-2.2 


3.06x102'' 


-3.49 


-1.3 


2.57x1030 


Ni 


0.00 


0.0 


2.16x10" 


0.13 


0.1 


3.52x1031 


Cu 


-3.06 


-1.1 


4.68x103° 


-2.08 


-0.7 


2.46x1028 


Zn 


0.17 


0.1 


2.93x10" 


-1.49 


-0.5 


6.46x102' 


Ga 


-2.09 


-0.7 


3.09xl02* 


-0.83 


-0.3 


3.54x1025 


Ge 


-1.48 


-0.5 


6.31x102'' 


-0.47 


-0.1 


3.11x1026 



the thermal expansion time scale. The different expansion 
timescales also result in a compositional asymmetry in ve- 
locity space. 

It is difficult to establish the observable features of our 
model since we have not conducted radiative transfer calcu- 
lations to generate light curves and spectra. It is currently 
unclear how much of an observational effect this asymme- 
try will have. A series of synthetic spectra generated over a 
range of time allows for direct comparison with observed su- 
pemovae. We conjecture, that even if compositional effects 
are obscured the spectra will show some dependence on ob- 



serving angle. This is because the side of the remnant that 
expands at higher velocities will also be at a lower density, 
making it more transparent at earlier times. 

Even though it is difficult to determine observational fea- 
tures of our model, it is instructive to compare and contrast 
our model with other recent results. The 'toy model' in 
iHille brandt et al l ( 2007) and Sim etjD (12007 !) is constructed 
similarly to our calculated yield distribution. They found that 
off-center distributions of burned material are likely to leave 
detectable imprints on observed light curves. An angular de- 
pendence of the light-curve peak brightness is introduced that 
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Fig. 6. — Final mass fraction of elemental(nonradioactive) Ni as a function 
of the ejection angle relative to the center of the star. The detonation started 
in the suiface layer of the star in the theta = 90° direction. The particles with 
a Ni mass fraction above 10^' are particles that have burned to NSE. 
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^ 10000 km s^^ faster on the side of the star opposite of the 
detonation. The ejecta in our model, however, do not have 
a sharp transition in velocity but a gradual increase in veloc- 
ity starting at a point ~ 90° away from the ignition point. It is 
difficult to determine if our model can reproduce the observed 
velocity gradients. It is worth noting that the deflagratio n to 
detonation transition model used in jMaeda et al.ll20ldl) did 
not produce the observed velocity shift or 'opening angle' ei- 
ther 

iKasen & Plewal ( 120051) attempted to calculate the spectral 
signatures of GCD by considering ejecta interacting with an 
extended atmosphere. The ejecta in their model were calcu- 
lated from a 1 -D simulation, and therefore lack the asymme- 
try in the nucleosynthesis that the react ion network calcula- 
tions of our 2-D simulation find. In ( Kasen & Plewal [20071) 
nucleosynthesis was done approximately with a 13 element 
reaction network. The surface flow, which consists partly of 
deflagration ash, which was excluded from our present model, 
needs to be considered. The surface flow might also have 
a spectral signature itself such as the presence of an high- 
velocity calcium absorber (iKasen & Plewall2005l) and should 
be compared with the underlying compositional asymmetry. 
In the case of sub-Chandrasekhar models, detonation of a 
pure helium shell leads to a layer containing iron-group el- 
ements such as titanium and chromium around the cor e ejecta 
(lFinketalJl20ia iSim et alllloTol: iKromer et aUlMoh . 

Another item to consider is that we do not include the ef- 
fects of metallicity on the nucleosynthesis. Our model is ini- 
tially composed of '^C and '^O. However, it has been shown 
that prior to the explosion of a carbon-oxygen white dwarf in 
a SN la there is a long period during which some '^C is con- 
verted into ^■'C as well as heavier elements (IChamulak et alj 
l2008h . This process makes even the most metal poor SNe la 
have a composition of more diverse than pure '^C and '^O. 
It is worth mentioning that for deflagration-detonation tran- 
sition (DDT) models where the detonation density was al- 
lowed to vary in relation to the flame speed as a function of 
metallicity dChamulak et al...2007,) the vield of ^^Ni produced 
also v aried with metallicity dBravo et"ani2010l: iJackson et al.l 
|20TO|) . DDT models wifli' varying metallicity and fixed det- 
onation density, howe ver, showed little variation in outcome 
jTownslev et al.ll2009h . 



Fig . 7 . — Final mass fraction of elemental iron group elements as a function 
of the final radial velocity. The particles with a Ni mass fraction above 10^^ 
are particles that have bumed to NSE. Particles ejected between —80° < 9 < 
—70° reach the highest radial velocities. 

might explain some over luminous SNe la. 

Our model has no deflagration ash so comparisons with de- 
flagration to detonation m odels, like [ Liv ne (1999), is prob- 
lematic. We can compare to' Maeda et al.l(i201 0). who showed 
that expansion velocity gradients as inferred from the Si II 
A 6355 absorption feature could be explained by a velocity 
shift of 3500 kms^' in Si happening in an 'opening angle' 
105 — 110° away from the ignition points caused by defla- 
gration ash. Even though our model has no deflagration ash 
a similar effect occurs with ^'^Si having ejection velocities 



The authors wish to acknowledge that this work was sup- 
ported by the US Department of Energy, Office of Nuclear 
Physics, under contract DE-AC02-06CH1 1357. This work 
was also supported by the NSF, grant AST-0507456, by 
the Joint Institute for Nuclear Astrophysics at MSU under 
NSF-PFC grant PHY 0822648, by the DOE through grant 
08ER41570, and by the Deutsche Forschungsgemeinschaft 
via the Emmy Noether Program (RO 3676/1-1). The soft- 
ware used in this work was in part developed by the DOE- 
supported ASC / Alliance Center for Astrophysical Ther- 
monuclear Flashes at the University of Chicago. Simula- 
tions presented in this work were run at the High Performance 
Computing Center at Michigan State University. 



REFERENCES 



Bravo, E., Dommguez, I., Badenes, C, Piersanti, L., & Straniero, O. 2010, 
ApJ, 711,L66 



Calder, A. C, Townsley, D. M., Seitenzahl, 1. R., Peng, F, Messer, O. E. B., 
Vladimirova, N., Brown, E. F, Truran, J. W., & Lamb, D. Q. 2007, ApJ, 
656, 313 



8 



Chamulak et al. 



Caughlan, G. R. & Fowler, W. A. 1988, At. Data Nucl. Data Tables, 40, 283 
Chamulak, D. A., Brown, E. R, & Timmes, F. X. 2007, ApJ, 655, L93 
Chamulak, D. A., Brown, E. R, Timmes, F. X., & Dupczak, K. 2008, ApJ, 
677, 160 

Cyburt, R. H., Amthor, A. M., Perguson, R., Meisel, Z., Smith, K., Warren, 
S., Heger, A., Hoffman, R. D., Rauscher, T., Sakharuk, A., Schatz, H., 
Thielemann, R K., & Wiescher, M. 2010, ApJS, 189, 240 

Fink, M., Ropke, R K., Hillebrandt, W., Seitenzahl, I. R., Sun, S. A., & 
Kromer, M. 2010, A&A, 514, A53+ 

Fryxell, B., Olson, K., Ricker, P., Tinunes, F. X., Zingale, M., Lamb, D. Q., 
MacNeice, P., Rosner, R., Truran, J. W., & Tufo, H. 2000, ApJS, 131, 273 

Fuller, G. M., Fowler, W. A., & Newman, M. J. 1982, ApJS, 48, 279 

Graboske, H. C, Dewitt, H. E., Grossman, A. S., & Cooper, M. S. 1973, 
ApJ, 181, 457 

Hillebrandt, W. & Niemeyer, J. C. 2000, ARA&A, 38, 191 
Hillebrandt, W., Sim, S. A., & Ropke, F. K. 2007, A&A, 465, L17 
lUadis, C, D'Auria, J. M., Starrfield, S., Thompson, W. J., & Wiescher, M. 

2001, ApJS, 134, 151 
Jackson, A. R, Calder, A. C, Townsley, D. M., Chamulak, D. A., Brown, 

E. R, & Timmes, R X. 2010, ApJ, 720, 99 
Kasen, D. & Plewa, T. 2005, ApJ, 622, L41 
— . 2007, ApJ, 662, 459 

Khokhlov, A. M. 1983, Soviet Astronomy Letters, 9, 160 

— . 1991, A&A, 245, 114 

Kromer, M., Sim, S. A., Fink, M., Ropke, F K., Seitenzahl, I. R., & 

Hillebrandt, W. 2010, ApJ, 719, 1067 
Langanke, K. & Martlnez-Pinedo, G. 2001, At. Data Nucl. Data Tables, 79, 

1 



Livne, E. 1999, ApJ, 527, L97 

Livne, E., Asida, S. M., & Hoflich, R 2005, ApJ, 632, 443 
Maeda, K., Benetti, S., Stritzinger, M., Ropke, F. K., Folatelli, G., 

SoUerman, J., Taubenberger, S., Nomoto, K., Leloudas, G., Hamuy, M., 
Tanaka, M., Mazzali, R A., & Elias-Rosa, N. 2010, Nature, 466, 82 
MazzaH, P. A., Ropke, F. K., Benetti, S., & Hillebrandt, W. 2007, Science, 
315, 825 

Meakin, C. A., Seitenzahl, I., Townsley, D., Jordan, G. C, Truran, J., & 

Lamb, D. 2009, ApJ, 693, 1188 
Plewa, T., Calder, A. C, & Lamb, D. Q. 2004, ApJ, 612, L37 
PodsiadlowsM, P., MazzaU, P., Lesaffre, P., Han, Z., & Forster, F. 2008, New 

A Rev., 52, 381 

Ropke, R K., Woosley, S. E., & Hillebrandt, W. 2007, ApJ, 660, 1344 
Seitenzahl, 1. R., Meakin, C. A., Lamb, D. Q., & Truran, J. W. 2009a, ApJ, 
700, 642 

Seitenzahl, I. R., Meakin, C. A., Townsley, D. M., Lamb, D. Q., & Truran, 

J. W. 2009b, ApJ, 696,515 
Sim, S. A., Ropke, F. K., Hillebrandt, W., Kromer, M., Pakmor, R., Fink, 

M., Ruiter, A. J., & Seitenzahl, I. R. 2010, ApJ, 714, L52 
Sim, S. A., Sauer, D. N., Ropke, R K., & Hillebrandt, W. 2007, MNRAS, 

378,2 

Townsley, D. M., Jackson, A. P., Calder, A. C, Chamulak, D. A., Brown, 

E. R, & Timmes, R X. 2009, ApJ, 701, 1582 
Woosley, S. E., Amett, W. D., & Clayton, D. D. 1973, ApJS, 26, 231 
Woosley, S. E. & Weaver, T. A. 1994, ApJ, 423, 371 



